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Abstract 

We construct an anomalous dimension for small x evolution which goes beyond standard 
fixed order perturbative evolution by including resummed small x logarithms deduced from 
the leading order BFKL equation with running coupling. Surprisingly, we find that once 
running coupling effects are properly taken into account, the leading approximation is very 
close to standard perturbative evolution in the range of x accessible at HERA, in overall 
agreement with the data, with no need for phenomenological parameters to summarise 
subleading effects. We also show that further corrections due to subleading small x log- 
arithms derived from the Fadin-Lipatov kernel can be kept under control, but that they 
involve substantial resummation ambiguities which limit their practical usefulness. 
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1. Introduction 

In recent years the theory of scahng violations for deep inelastic structure functions at small 
X has attracted considerable interest, prompted by the experimental information coming 
from HERA [1,2]. New effects beyond the low-order perturbative approximation [3-5] to 
anomalous dimensions or splitting functions should become important at small-a;. The 
BFKL approach [6-10] provides in principle a tool for the determination of the small- 
X improvements of the anomalous dimensions [11,12]. However, no major deviation of 
the data from a standard next-to-leading order perturbative treatment of scaling viola- 
tions [13,14,1,2] has been found, and a straightforward [15] inclusion of small-x logarithms 
is in fact ruled out by the data [16,17]. By now the origin of this situation has been mostly 
understood [18-21]. 

The BFKL kernel x{^s, M), as is well known, has been computed to next-to- leading 
accuracy (NLO): 

x{M,as) = asXo{M) + alxi{M) + .... (1.1) 

The problem is how to use the information contained in xo and xi in order to improve 
the splitting function derived from the perturbative leading singlet anomalous dimension 
function ^{as, N) which is known [4] up to NLO in ag'- 

7(iV,a,)=a,7o(iV) + ahi(^) + (1-2) 

in such a way that the improved splitting function remains a good approximation down to 
small values of x. This can be accomplished by exploiting the fact [22] that the solutions 
of the BFKL and GLAP equations coincide at leading twist if their respective evolution 
kernels are related by a "duality" relation. In the fixed coupling limit the duality relation 
is simply given by: 

xiliN,as),as) = N. (1.3) 

The splitting function then will contain all relative corrections of order (ct^ log 1/x)", de- 
rived from xo{M), and of order Q!s(q;s log 1 /a;)", derived from xi(M). 

The early wisdom on how to implement the information from xo was completely 
shaken by the much softer behaviour of the data at small x and by the computation of 
Xi [8-10], which showed that the naive expansion for the improved anomalous dimension 
had a hopelessly bad behaviour. In refs. [20,19] we have shown that, as a consequence 
of the physical requirement of momentum conservation, these problems are cured if the 
small-a; resummation is combined with the standard resummation of coUinear singularities, 
by constructing a 'double-leading' perturbative expansion. However, the next-to-leading 
correction remains large, and it qualitatively changes the asymptotic small-x behaviour of 
structure functions by changing it from x~^° to x^'^ — x~^"e^^^ ~ x~^° [1 -|- AA^ -|- ....], 
with: 

A = Ao + AA, Xo = asxoii) = oisco, AA = a^xi(|) + • • • • (1-4) 

The effects of this correction can be resummed, but there are ambiguities in the procedure 
which entail the dependence of the results on free parameters, or, alternatively, model 
assumptions. 
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A traditionally delicate point of the small-a; resummation is the treatment of the run- 
ning of the coupling [23-28]. While it has been known for some time [29] that running 
coupling effects can be included perturbatively order by order at small x, in a recent pa- 
per [30] we have shown that an all-order resummation of running coupling effects at small 
X is desirable and can in fact be accomplished. Indeed, in a perturbative approach, running 
coupling terms can be included by adding effective subleading A^j contributions to the 
BFKL kernels xi, X2-, ■ ■ ■, which turn out to have singularities at M = 1/2. These singular- 
ities correspond [29,30] to an enhancement by powers of Inl/a; of the associated splitting 
functions, which may offset the perturbative suppression by powers of ag. In ref. [30] we 
have shown that the all-order resummation of these contributions is fully compatible (up 
to higher twist corrections) with the standard factorized perturbative evolution of parton 
distributions, and can thus be performed at the level of splitting functions. The resummed 
splitting functions remain smooth in the small-x limit, despite the well-known [23-26] fact 
that the naive running coupling BFKL solution instead displays oscillatory instabilities at 
small X. 

The purpose of this paper is to discuss the physical impact of this running coupling 
resummation, and describe in detail its phenomenological implementation. We will argue 
that we now know the way the information contained in xo{M) should be used in order to 
construct a better first approximation for the improved anomalous dimension. Indeed, we 
find that, once running coupling effects are properly included in the improved anomalous 
dimension, the asymptotic behavior near x = is much softened with respect to the naive 
Lipatov exponent. Hence, the corresponding dramatic rise of structure functions at small 
X, which is phenomenologically ruled out, is replaced by a milder rise, whose steepness is 
determined by the Lipatov exponent. It then follows that a leading-order approximation 
based on the standard BFKL kernel xo is phenomenologically viable. 

Our proposed leading approximation is of a reasonably simple explicit form, directly 
suitable for the determination of scaling violations and fits to the data, contains no free 
parameters, and is in agreement with the general trend of the data. We find it remarkable 
that such a first approximation indeed exists and apparently works so well, so that we 
devote a good fraction of this article to describing it and its physical foundation. We then 
construct and discuss a perturbative expansion based upon it, and show in particular that 
the next order correction, which includes the effects of the subleading kernel Xi, leads to 
small corrections for a reasonable range of the resummation parameters. However, some 
parameter dependence and resummation ambiguities are necessarily present at this level, 
and in fact at next-to-leading order the associated ambiguity is of the same size as the 
correction itself. The fact that the corrections are small but ambiguous justifies the use of 
the simple leading-order approximation for practical purposes. 

2. Notation and Basic Formalism 

The behaviour of structure functions at small x is dominated by the large eigenvalue 
of evolution in the singlet sector. Thus we consider the singlet parton density 

^(e, t) = x[gix, Q^) + k,® qix, Q% (2.1) 
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where ^ = logl/x, t = logQ'^/jj^, g{x,Q'^) and q{x,Q'^) are the gluon and singlet quark 
parton densities, respectively, and kq is such that, for each moment 

G{N,t)= / (i^e-^^G(e,t), (2.2) 
Jo 

the associated anomalous dimension '^{as{t),N) corresponds to the largest eigenvalue in 
the singlet sector. The generalization to the full two-by-two matrix of anomalous dimen- 
sions is discussed in detail in ref. [20]. 

At large t and fixed ^ the evolution equation in A'"-moment space is then 

jG{N, t) = 7(«.(t), N)GiN, t), (2.3) 

where as{t) is the running coupling. The anomalous dimension is completely known at 
one- and two-loop level, as given in eq. (1.2). The corresponding splitting function is 
related by a Mellin transform to ^{as, N): 

^{as,N)= f dxx^P{as,x). (2.4) 
Jo 

Small X for the splitting function corresponds to small for the anomalous dimension: 
more precisely P ~ l/a;(log(l/x))'^ corresponds to 7 ~ n\/N'^^^ . Even assuming that a 
leading twist description of scaling violations is still valid in some range of small x^ as 
soon as x is small enough that ctsC ~ I5 with ^ = logl/x, all terms of order {as/x){asC)^ 
(LLx) and (y.s{(y.s/x){as^)'^ (NLLx) which are present in the splitting functions must be 
considered in order to achieve an accuracy up to terms of order a1{as/ x){as^)^ (NNLLx). 

As is well known, these terms can be derived from the knowledge of the kernel 
x(cts, given in NLO approximation in eq. (1.1), of the BFKL ^-evolution equation 

^G(^,M) = x(««,M)G(e,M), (2.5) 

which is satisfied at large ^ by the inverse Mellin transform of the parton distribution 

/oo 
dte-^*G(e,t). (2.6) 
-00 

This derivation was originally performed [11] at LLx by assuming the common validity of 
eq. (2.5) and eq. (2.3) in the region where and ^ are both large. However, it was more 
recently realized [22,19,20] that the solution of eq. (2.5) coincides generally with that of 
eq. (2.3), up to higher twist corrections, provided only that the kernel of the former is 
related to that of the latter by a 'duality' relation, and boundary conditions are suitably 
matched. This implies that the domains of validity of these two equations are in fact the 
same in perturbation theory, up to power-suppressed corrections. 
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The derivation of duality is simplest when the coupling does not run, in which case 
the relation between the kernels of the two equations is given by eq. (1.3). A simple proof 
of this statement was given in ref. [19]. Expanding ^(a.si N) in powers of cug at fixed ag/N 

7(a„ N) = -fsias/N) + a,7,,(a,/iV) + . . . , (2.7) 

and x(cKs) M) in powers of ag at fixed M 

X{as, M) = a,xo(M) + alxi{M) + . . . , (2.8) 

we then find that xo determines jsio^s/N), while agXi leads to ag^ssio^s/N): 

XoM^)) = -, (2.9) 



a 



s 



ry ( - _^^^S^IAjl}1 (2 10) 



Upon Mellin inversion, 'fs{as/N) and as'^ss{cts/N) correspond respectively to all terms of 
order {as/x){asC)^ and as{as/x){asO^ the splitting functions. 

When one goes beyond LLx, i.e. beyond the leading-order approximation for x, the 
running of the coupling cannot be neglected, and the duality relation must be re-examined. 
Indeed, in M space the usual running coupling as{t) becomes a differential operator and, 
by taking only the one-loop beta function into account, one has: 

fy 

S. = -— ^ + ---, (2.11) 

where (3o is the first coefficient of the /3-function (so (3 = —(3oa'^ H — ■ and (3o = 0.663146... 
for 4 fiavours), with the obvious generalization to higher orders. Hence, assuming the 
coupling to run in the usual way with Q^, the ^-evolution equation eq. (2.5) becomes [31] 

^G(e, M) = x(S., M)G(e, M), (2.12) 

where the derivative with respect to M acts on everything to the right, and x ™ay be 
expanded as in eq. (1.1) keeping the powers of on the left. In ref. [30] we have shown 
that eq. (2.12) can indeed be viewed as an alternative representation of the standard 
renormalization-group equation. 

It is clear from eq. (2.12) that running coupling effects begin at NLLx. To NLLx one 
finds [29] that the solution of eq. (2.12) is again the same as that of a dual t-evolution 
equation (2.3), provided the fixed-coupling duality relation eq. (1.3) is modified by letting 
as Ois{t), and then by adding to jss eq. (2.10) an extra term [23] A^ss proportional to 

^^ss{^) = -Po^^^. (2.13) 
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Equivalently, the duahty relation eq. (1.3) can be formally preserved, provided that cug 
ctsit) and the function x used in it is no longer identified with the BFKL kernel, but rather 
given by an 'effective' x function 

Xeff(a.,M) = x(as,M) + Ax(a.,M), 

(2.14) 

Ax{as,M)^alAxi{M) + ..., 

where xi^sj M) is obtained letting — > ag in the kernel of eq. (2.12), and the correction 
term Ax to NLLx is given by 

The problem with this perturbative approach is that the correction terms which must 
be included in the effective x (2-14) have an unphysical singularity: xo{M) has a 
minimum at M = so the denominator of Ax;i(M) vanishes, resulting in a simple pole 
in the NLLx correction Axi eq. (2.15) at M = |. The NNLLx correction Ax2 turns out 
to have a fourth-order pole, and in fact at each extra order three extra powers of (xo)""^ 
appear. As explained in ref. [29], this leads to a perturbative instability: as a consequence 
of the singularity, the splitting function AP^s eq. (2.4) associated with the anomalous 
dimension Ajgg behaves as 

p . t^ ~ («sO • (2-16) 

This problem was completely solved in ref. [30] , where we have shown that the M = | 
singularity and the corresponding ones that appear at higher orders in agPo are an ar- 
tifact of the expansion, are not present in the all-order solution and can thus be elimi- 
nated, thanks to the possibility of matching the all-order solution to the perturbative one. 
The corresponding resummation, which will be discussed in detail in the next section, is 
obtained by explicitly solving eq. (2.12) for xi^s^M) ~ cisXo, with xo replaced by its 
quadratic approximation near M = |. As is well known, the solution is given in terms of 
an Airy function. This solution is added to the anomalous dimension with subtraction of 
the NNLx expansion terms to avoid double counting. The x function dual to the resulting 
expression for 7 is now regular near M = |. 

In the sequel, we will discuss in detail the form and the properties of the improved 
anomalous dimension. In the next section we describe the leading approximation, where as 
input we only use the one- loop perturbative anomalous dimension as'^oi^) the leading 
order BFKL kernel asXo{M) but with the running coupling effects taken into account. 
Subleading corrections, including those from xi, will then be discussed in section 4. 



3. The Improved Anomalous Dimension: First Iteration 

Assuming that one only knows 70 (AT) and xo{M), we argue that the optimal use that 
one can make of these inputs is to write down for the improved anomalous dimension the 
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following expression: 



7/(a„ N) = K7o(iV) + 7,(^) - ^] + 

(3.1) 



+ 



7a(co, as, iV) - I + ^ ^^^[N - agCo] + jPoas\ - mom. sub. 



The first line on the right-hand side, within square brackets, is the usual double-leading [19] 
expression for the improved anomalous dimension at this level of accuracy, made up of the 
first order perturbative term q;s7o plus the power series of terms (as/N)"' contained in 
7s(^), obtained from xo using eq. (2.9), with subtraction of the order as term to avoid 
double counting {ca = Uc = 3). In the second line, the second pair of square brackets 
contain the running coupling Airy term: 7^(00,0:5,-^) is the Airy anomalous dimension, 
and the remaining terms subtract the contributions to 7^(00, cts,-^) which are already 
contained in 7^ and 70 as we shall explain shortly. Finally "mom. sub." is a subleading 
subtraction that ensures momentum conservation ji{as,N — 1) = 0. We now specify 
7a(co, CKs, AT) and "mom. sub.". 

For 7a(c, as, N) ^ we start from the Mellin transform of eq. (2.12), namely [30] 

NG{N, M) = x{as, M)G{N, M) + F(M), (3.2) 

with kernel 

Xq{as, M) = as[c+ ^«(M - Msf]. (3.3) 

The expression within squared brackets is the quadratic approximation of x near its min- 
imum M — Ms- At leading order, Mg, c (the Lipatov exponent) and k are given by: 

M, = M« = i, (3.4) 

477, 

c = Co = —log 2 = 2.64763..., (3.5) 

TT 

K = Ko = -2^V'"(|) = 32.1406..., (3.6) 

where il){x) is the digamma function. The Airy anomalous dimension is found from the 
(Mellin transformed) solution Gq{N,t) of this equation: 

..(c, N) . 4 lnG,(iV, «) = M. + ('-^\ 4S^t^. (3-7) 



dt ' ' ' \ K J Ai[z{as{t),N)] ' 

where Ai{z) is the Airy function, which satisfies 

Ai"{z) - zAi{z) = 0, (3.8) 



^ For convenience we only include c (and not k, see below) among the explicit arguments of 
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with Ai(0) = 3-2/3/r(2/3), 



z{as{t),N) 



2f3oNY^^ 1 



C 



as{t) N 



(3.9) 



Along the positive real axis, the Airy function is a positive definite, monotonically 
decreasing function of its argument, and it behaves asymptotically as 

Ai(^) = ^n-'/h-'/''exp{-lz^/^){l + 0{z-^/^)). (3.10) 

Thus the large-2; behaviour of 7^, which corresponds to N large and/or as{t) small, is 
given by 



7a(c, as, N) = if{as/N) + ^ ^ 



(3.11) 



\ K las 

The leading 7^ behaviour coincides with the 'naive', i.e. fixed-coupling dual anomalous 
dimension: namely, the anomalous dimension which is found from Xq sq. (3.3) using the 
leading-order duality relation (1.3). F\irthermore, the leading 0{al/N) correction is sim- 
ply obtained from eq. (3.10): 

as^fsMN) = ---^ (3.12) 
41- 

and it coincides with the next-to-leading running coupling correction eq. (2.13) computed 
for the quadratic kernel eq. (3.3). Subsequent terms in the asymptotic expansion eq. (3.11) 
correspond to the further a'^jsss, ■ ■ ■ corrections which are found when the perturbative 
solution of the running coupling eq. (2.12) is pursued to higher orders. The asymptotic 
series is (Borel) resummed [30] by the Airy anomalous dimension. At leading order, the 
double-counting between the double-leading terms and Airy anomalous dimension is re- 
moved by subtracting the leading term of the asymptotic expansion eq. (3.12), as well as 
the term q;s7^(0) = —^Poas, which is 0{as) and thus must be subtracted because it is 
already in asjo- 

The momentum conservation subtraction "mom. sub." is an 0{a^) term which can 
be defined as: 

mom. sub. = g{N)^i{as, 1), (3.13) 

where g{N) is a weight function with ^(1) = 1, with no singularities for A?" > and 
'yi{as,N) is given by 

7/(q;s, N) — 7/(q;s, N) — mom. sub.. (3-14) 



We can take g{N) = 1 or, for instance. 
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Figure 1: The improved anomalous dimension 'yi(as,N) for c = Co, k = kq, Ms = 1/2 
(i.e. the values in eqs. (3.5), (3.6), (3.4)that are obtained from xoiM)), g{N) = 1 and 
tts = 0.2 (solid), compared with 7dl-lo (the curve with a step) and GLAP LO and 
NLO (dashed and dotted). 



where r is a positive real number. The possible advantage of this choice is to have a 
damped large N behaviour. These different choices are all equivalent up to subleading 
terms. 

We now discuss the properties of the improved anomalous dimension in this approxi- 
mation. In the limit as ^ with arbitrary N, 7/(as, A^) reduces to ag'joiN) + 0{a^): in 
particular as — oo the Airy term reduces to an 0{a1) constant (the N oo limit of 
the next term, 7^^, in the expansion of eqs. (3.10)-(3.12) up to an 0{N~^) correction). 
For as ^ with as/N fixed, -fi{as,N) reduces to 7dl-lo = as7o(A^) + ls{j^) - 
i.e. the leading term of the double-leading expansion. Thus the Airy term is subleading in 
both limits. In spite of this, its role is very significant because of the singularity structure 
of the different terms in eq. (3.1). In fact, 70 (-^) has a pole at = 0, 7s has a branch 
cut at = CKsCo, and 7^ has a pole at A^ = A'o < «sCo, where A'o is the position of the 
rightmost zero of the Airy function. The importance of the Airy term is that the square 
root term subtracted from 7^ cancels, within the stated accuracy, the branch cut of 7^ at 
A^ = agCo and replaces the corresponding asymptotic behaviour at small x with the much 
softer one from 'Ja- 

The effect of this replacement is clearly seen from figure 1, where we compare the plot 
of 7/ (solid curve) as a function of A^ with those of 7dl-lo (the curve with a step, starting 
at the cut branch point) and of GLAP LO and NLO (which are nearly superimposed). The 
curves are computed with c = cq, k = kq, Ms = 1/2 (i.e. the values in eqs. (3.5)- (3.4) 
that are obtained from xoiM)), to the momentum subtraction with g{N) = 1 eq. (3.15) 
and CKs = 0.2, a value which is roughly appropriate for the HERA kinematic range. Since 
A'o is significantly smaller than ctgCo, the improved anomalous dimension 7/ remains close 
to the GLAP curves down to much smaller values of A^ than 7dl-lo- It is only at A^ ^ 0.3, 
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Figure 2: The improved splitting function Pgg{x) for c = Co, k = kq, Ms = 1/2 (i.e. the 
values in eqs. (3.5), (3.6), (3.4) that are obtained from xoiM)), g{N) = 1 and as = 0.2 
(solid), compared with that from 7dl-lo (dot-dashed) and GLAP LO (dashed). 



where 77 > 0.5, that the improved anomalous dimension becomes sizably different than 
the NLO perturbative anomalous dimension (for the chosen value of ag the pole of 7/ is 
at No = 0.211237.. .while a^co = 0.529525...). 

The deviation of 7/ from GLAP at small values of has a moderate effect on the 
splitting function. In fig. 2 we plot Pgg{x) corresponding to 7/ (dashed) compared to the 
perturbative LO (solid), while the rising curve at small x corresponds to 7dl-lo- This 
plot makes particularly clear that the improvement obtained by implementing the Airy 
asymptotics is really important at small x, especially when we take into account that the 
data seem to follow rather closely the unresummed perturbative evolution: the rapid small 
X rise of the splitting function is removed by the resummation, and in fact the resummed 
splitting function closely follows the perturbative one down to x ~ 10"'^. 

The improvement of the agreement with the leading order GLAP anomalous dimen- 
sion when the double-leading anomalous dimension is supplemented by running coupling 
resummation according to fig. 1 follows from three reasons. First, the location of the right- 
most singularity is lowered from ctgCo to A'o- Second, the double-leading branch cut is 
replaced by a pole, which is the same type of singularity as in GLAP. Finally, the residue 
of the pole in is more than an order of magnitude smaller than the residue of the GLAP 
pole. Indeed, a straightforward calculation using eq. (3.7) shows that 

7AK(t), N) = J^^JA + 0[{N - Nof] 

When as = 0.2, = 0.014... to be compared with the value 70 = 3as/njj + 0[N^] = 
0.191.../A^ + 0[A^'^] of the residue of the leading-order GLAP pole. The combination of 
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these facts implies that the singularity in 7^ only kicks in rather late, i. e. for rather small 
values of A^. This is apparent in fig. 2, where the rise rAX~^° due to the Airy pole is only 
visible at very small x on top of the constant behaviour related to the GLAP pole. Due 
to momentum conservation, which fixes the integral of xP, the rise must be compensated 
by a small depletion of the splitting function in the intermediate x range, also visible 
in the figure. While all these properties are generic, the quantitative agreement of the 
resummed result with leading order GLAP is only good when the intercept is close to 
the leading order value ctsCo [eq. (3.5)] determined from the BFKL kernel xo- Running 
coupling resummation is the essential ingredient that reconciles the BFKL intercept with 
the GLAP behaviour of the data GLAP [14,32]. 

In our previous work in refs. [19,20] we had adopted a more phenomenological approach 
in order to cope with the fact that the rise of the singlet structure function implied by the 
Lipatov exponent (corresponding to the anomalous dimension with the cut at ctgCo in fig-1) 
is too sharp in comparison to what is seen in the data. Noting that the correct asymptotic 
exponent is not reliably determined by the two known terms of its perturbative expansion 
(in fact the computed correction to the leading result is quite large), we had treated it 
as a parameter to be fitted from the data. To this effect, we interpreted the parameter 
A A in eq. (1.4), formally of order parameter to be fixed empirically. The lack of 

any indication of a power-rise in the data suggests then that AA should be negative and 
possibly numerically as large as, say, 0{as). In leading order, the corresponding improved 
anomalous dimension becomes 

lioldias, N) = [asio{N) + ls{j^) - ^] - mom. sub. (3.17) 

The shift (N — AA) in 7s provides the change from x~'^° to x~^ = a;~^°e^^^. As 

discussed in ref. [19], there is an ambiguity in the subtraction ncas/{nN) because there 
too we could make the replacement N ^ {N — AA). In eq. (3.17) we have adopted the 
subtraction with AT, denoted as R-subtr action in ref. [20], because it is more suitable for 
the present purpose of comparing the old approach with 77 in eq. (3.1), as we shall see 
shortly. 

What we have accomplished in this section is to show that the resummation through 
the Airy asymptotics of higher-order contributions to the anomalous dimension related to 
the running coupling corrections is by itself sufficient to produce the required softening 
of the behaviour at small x. Thus it is no longer necessary to introduce A A by hand 
as a free parameter (at least at leading accuracy). In fact, we had already found in 
ref. [19] that the value of AA that provides the best matching to GLAP is exactly such 
that A = CKgCo + AA ~ 0.21 (for as ~ 0.2) in remarkable agreement with the value 
A^o = 0.211237 determined here from the Airy asymptotics. The fact that the data are in 
good agreement with unresummed perturbative evolution suggests that the true all-order 
value of A is likely to be quite close to the leading-order value, so AA is presumably rather 
small and there is perhaps no need to introduce it as a free parameter in order to fit the 
data. 

The sensitivity of our new approach to the value at the minimum A = CKgC = agCo + AX 
can be demonstrated by going back to eq. (3.1), replacing N hy N — A A in 7s and cq by 
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Figure 3: The effect on the improved anomalous dimension of varying the value of x 
at the minimum: A = OsC = OsCo + AA, as described in the text, with as = 0.2 and 
ttaC = 0.4, 0.529525..., 0.6. The dashed curve is the LO GLAP curve, the solid curve 
corresponds to the BFKL value a^c = 0.529525... while the upper and lower (at small 
A'^) dot-dashed curves refer to otsC — 0.6, 0.4, respectively. 



c in the Airy term (including the square root subtraction): 



77(a., iV, c) = [a.7o(iV) + ls{j^) - ^] + 

I (3.18) 

+ 7a(c, Us, N)-\^ \J l^S^ ~ (^sc] + \(3Qas - mom. sub. 

Note that we have not performed the AA shift in the double-counting subtraction term, 
corresponding to the R-subtraction of ref. [20] . This is necessary, otherwise the cancellation 
of singularities at = agC between the double-leading anomalous dimension and the Airy 
term would be spoiled. The resulting improved anomalous dimension for = 0.2 and 
ctgC = 0.4, 0.529525... (the BFKL value, see eq. (3.5)), and 0.6 are shown in figure 3. The 
corresponding splitting functions are shown in figure 4. We see that changing c in this 
range does not alter the overall trend of the anomalous dimension. However the impact on 
the splitting function at small x is quite noticeable and, in fact, the observed approximate 
validity of the perturbative fits suggests that the viable range of AA is inside the interval 
presented in figure 4. This possibility of varying c in the leading term could be adopted 
as a phenomenological way to optimize the leading formula on the data. We have also 
checked that moderate variations of and Mg eqs. (3. 4), (3. 6), do not alter the general 
picture. 

Since the reference perturbative fits are based on the NLO GLAP it may be desirable 
to modify our improved anomalous dimension in such a way that at small ctg with ag/N 
fixed it reduces to NLO GLAP given by eq. (1.2): 7(A^, as) = as7o(A^) + a'iji{N). This 
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Figure 4: The effect on the improved splitting function of varying the value of x ^■t 
the minimum: A = OsC = UsCq + AA, as described in the text, with as = 0.2 and 
OsC = 0.4, 0.529525..., 0.6. The dashed curve is the LO GLAP curve, the solid curve 
corresponds to the BFKL value asC = 0.529525... (same as the solid curve of Fig. 1) 
while the upper and lower (at large 1/x) dotdashed curves refer to a^c = 0.6, 0.4, 
respectively. 



can readily be obtained by modifying eq. (3.1) in the following way: 

7r(«„ N) = K7o(iV) + al^^{N) + 7,(^) - ^] + 

+ 7^(co, tts, + [N - a.coj (3.19) 

(y. 

+ \(3oas{l + ^co) - mom. sub. 

Here we used the fact that in the expansion of 7s(as/A^) in powers of its argument the 
quadratic term is absent, so that it is sufficient to only add a subtraction to the Airy term 
in order to guarantee that in the perturbative limit the NLO GLAP anomalous dimension 
is exactly recovered. This subtraction (the first term in the third line of eq. (3.19)) is 
obtained by expanding as7^ eq. (3.12) to order a1 at fixed A^. This formula improves 
the NLO GLAP with the leading terms of the BFKL corrections with running coupling. 
The result is shown in figure 5 and compared to the GLAP anomalous dimension and 
the resummed LO result of Figure 1. It appears that, thanks to the extra next-to- leading 
subtraction in eq. (3.19), the effect of the resummation is further reduced, in particular 
correcting the slight distortion of the LO resummed result in the 0.5 ^ N ^ 1.5 region 
seen in the figure (dotdashed curve). 

In conclusion we have shown that, given the tendency of the data to closely follow 
the GLAP evolution down to the smallest values of accessible to HERA, the most 
effective way to implement the information contained in the leading BFKL kernel xo is to 
use duality, the double-leading prescription (that improves xo near M = in agreement 
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Figure 5: The improved anomalous dimension -^f" eq. (3. 19) (solid), LO (dotted) and 
NLO GLAP (dashed), which is the corresponding perturbative limit in this case, for 
c = Co, K = Ko, Ms = 1/2 (i.e. the values in eqs. (3.5), (3.6), (3.4) that are obtained 
from xo(Af)), g{N) = 1 and as = 0.2. The resummed LO curve (solid curve of Figure 
1) is also shown for comparison (dotdashed). 

with momentum conservation), and the Airy asymptotics to get just the right amount of 
softening of the smaU x exponent. In the next section we wiU show that the non leading 
corrections, including those arising from xi, do not spoil this overall picture. 

4. Perturbative expansion of the resummed anomalous dimension 

In the previous section we have seen that at leading order the resummed anomalous 
dimension eq. (3.1) and the GLAP anomalous dimension are in very good agreement. 
However, this result is phenomenologically significant only if the resummed anomalous di- 
mension eq. (3.1) is stable upon higher order corrections. To investigate whether this is the 
case we must discuss how eq. (3.1) can be viewed as the leading order of a resummed per- 
turbative expansion of the anomalous dimension. As suggested in ref. [30] , this can be done 
by using the resummation of running coupling effects by means of the Airy asymptotics 
to improve the expansion eq. (2.7) of the anomalous dimension in powers of as at fixed 
ag/N, and then combining this improved expansion with the standard GLAP anomalous 
dimension, to obtain a running-coupling improvement of the double leading expansion, 
thereby generalizing to all orders eq. (3.1). 

Hence, we must consider the running-coupling resummation of the small-x expansion 
of the anomalous dimension eq. (2.7) beyond leading order. Higher order corrections Xi 
to the kernel xo will change in general the values of all the parameters (3.4)-(3.6) which 
characterize the quadratic approximation to x near its minimum. Whenever the value 
of X ci-t its minimum differs from the leading-order value, the perturbative expansion of 
the kernel x must be reorganized if we wish to obtain a stable expansion of the splitting 
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function, i.e. such that the ratios Pss/Pai Pass/Pai ••■ cill remain finite as ^ ^ oo. A 
necessary and sufficient condition for this to be the case is [29,19] that to each order a 
constant is subtracted from Xi a-nd added back into Xo^ iii such a way that the subtracted 
Xi no longer lead to a shift in the minimum of xo- At leading order, this amounts to 
performing the shift N ^ N — A A as in eq. (3.18). 

In view of a running coupling resummation, however, it is convenient to include in 
Xo the full leading asymptotic behaviour, by subtracting not only a constant, but also 
the leading linear and quadratic terms in the Taylor expansion of x about its all-order 
minimum. Namely, assume that the all-order x has a minimum in the vicinity of which it 
can be approximated by the quadratic expression of Xg, as given in eq. (3.3), where now 

M, = MO + a,MW + ..., (4.1) 
c = Co + asCi + (4.2) 
K = KO + CKsKi + (4.3) 

We then reorganize the expansion of x eq. (1-1) as follows: 
X(M, as) = asXo + «^xf + • • • ; 

xf = X^-xf, ^>0; 

X« = Xi{Ms) + (M - Ms) XiiMs) + |(M - Ms)'x'I{Ms), i > 0. 



(4.4) 



To any finite order, the expansion eq. (4.4) of x diff'ers from the standard one eq. (1.1) by 
subleading terms: 



n 

J2<+'Xim = $^«l+^xg^(M) + 0[<+i]. (4.5) 



i=0 i=0 



Therefore, eq. (4.4) is a bona fide resummed version of the original expansion. We may then 
use the resummed expansion to determine order by order 7^, 7^ and so forth, through 
perturbative duality with the suitable running coupling corrections. In particular, the 
next-to-leading order running coupling contribution to 7^ is still given by eq. (2.13), but 
with Xo replaced by its resummed counterpart Xo^ eq. (4.4). This leads to a resummed 
expansion of 7 in powers of as at fixed ag/N which, to any finite order n, differs by the 
standard expansion eq. (2.7) by subleading 0{a'^~^^) terms. Because, by construction, x^r 

do not shift the minimum of Xr^ 1 the ensuing resummed expansion of 7 is stable, just like 
the expansion discussed in our previous work refs. [29,19,20]. 

This new reorganized expansion differs from that discussed in refs. [29,19] because 
now not only the value of x at the minimum, but also its curvature is reshuffled into the 
leading order, and furthermore, the expansion is performed about the all-order minimum 
Ms instead of the leading order minimum M = |. The leading order resummed Xo^ has 
in general not only a value at the minimum c 7^ cq, like the leading order resummed xo 
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of refs. [29,19] but also a minimum at Mg ^ | and curvature k ^ kq. Only if k = kq 
and Mg = i, then Xo^ would be of the same form as the leading order resummed Xo of 
refs. [29,19]. Note, however, that even so, the subsequent terms xf" with i > 1 in the 
expansion of x would not be the same as the corresponding Xi of refs. [29,19] (compare 
eq. (4.4)). 

The resummed 7 at any desired order can now be combined with the Airy anomalous 
dimension 7^ determined from the quadratic kernel with generalized (4.1)-(4.3) values of 
the parameters. However, if the location of the minimum Mg differs from its leading order 
value, Mg ^ i, already at next-to-leading order this is no longer sufficient to resum all 
singularities related to the running coupling term A7SS eq. (2.13). Indeed, if Mg ^ |, Xo 
is no longer symmetric about M = Mg and therefore its Taylor expansion about M = Mg 
also contains in general a cubic term. It is straightforward to see that in such case A^gg 
acquires a new singularity. 

Consider the case of a generic kernel x, with a minimum at M = Mg, and containing 
terms up to fourth order: 



x(^)(a.,M) = xg + a 



l^^iM - Mgf + ^{M - MgY 



(4.6) 



Expanding the anomalous dimension 'fi'^\ag/N) which is dual (eq. (1.3)) to this x*^^"* about 
its rightmost singularity ag/N = c yields 




e 

k4! 



K \ag 



k3! V k \as 
+ 0[S']+0[e'] }, 



(4.7) 



so higher order terms in the expansion of x in powers of M — Mg correspond to higher order 
contributions in the expansion of its dual anomalous dimension in powers of 0v7 ag — c. 

Substituting 7^^"* (eq. (4.7)) in the expression eq. (2.13) of the running coupling contribution 
to 7ss leads to 



1 88 — 



x[)-(7l^^)Xo(7i^^) 



7f.+ 



M [2 



24 



K 



3/2 



N 




a. 



+ 



TV 



OLg 



(4.8) 



where 7^ is the running coupling contribution eq. (3.12) which is resummed by the Airy 
anomalous dimension. The cubic term induces the square-root singularity proportional to 
8 in eq. (4.8), which must also be resummed. Because the singularity is linear in 5, this 
resummation can be performed by treating the cubic term as a leading-order perturbation. 
Note that the quartic term does not induce any further singularities to this order. 

In order to resum running-coupling singularities which are associated to the cubic 
term in the expansion of x about its minimum, we determine the (Mellin transformed) 
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solution G^^\N,t) to eq. (2.12) when the cubic term is also included in the kernel: 

«-(-,0=(H^)"exp^e.p-{A(H^)'-^},,,.,,,.,. 

(4.9) 

Using the identity Ai"{z) = zAi{z) and expanding in powers of 5 we get the anomalous 
dimension 



7i'^(a.(t),iV) = |lnGf (iV,t) 



2 



5/3 



2z{a,{t),N)- 



Ai'{z{as{t),N)) 
Ai{z{as{t),N)) 



+ 0[d\ 
(4.10) 

where 7^ is given by eq. (3.7). 

The last term inside the square bracket in eq. (4.10) has a double pole at N = Nq, 
where Ai{z{as, N)) vanishes. This is due to the fact that the location of the zero of 
G'^^\N,t) eq. (4.9) is shifted away from A'^o by the 6 correction, and therefore the pole in 
the anomalous dimension is also shifted. To first order in 6, the zero moves from A^o to 

iVo = ^^0 + ^TA, (4.11) 

with rA given by eq. (3.16). Because ta is 0(as), the ensuing correction to the asymptotic 
behaviour is subleading [30].^ It follows that the linearized expression eq. (4.10) of the S 
correction to the anomalous dimension has a double pole in Nq, and is correct provided 
if N — Nq is not too small, i.e. if x is not too small. For the typical values of the 
parameters {agC ~ 0.4), Nq — Nq ~ 10~^, so the correction to the asymptotic behaviour 
x~^o _^ x~^o is approximately ~ 10% when x ~ 10~^, and only becomes of order 
one when x is extremely small. Hence, even though the linearized anomalous dimension 
eq. (4.10) does not reproduce the correct leading singularity of the full anomalous dimension 
determined from the solution eq. (4.9), it is a good approximation to it for all practical 
purposes. 

We can now determine the expansion of 7)^ in powers of as at fixed ag/N 

if = 7E,(«./iV) + c^slZsi^^/N) + (4.12) 
by using the asymptotic expansion eq. (3.10): 



lfM.IN)=-,}-l.(l\ f^-c], (4.13) 



12 \K J \a 



•^8 



^ Note that the explicit expressions of the correction to the solution and the location of the 
pole in the anomalous dimension given in eq. (5.2) and (5.3) of ref. [30] are incorrect, and should 
be respectively replaced by eq. (4.9) and eq. (4.11). The main conclusion that the modification 
of the asymptotic behaviour due to cubic and higher corrections is subleading remains however 
unaffected. 
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where 7^ and 7^ are given by eq. (3.11) and eq. (3.12) respectively. This verifies exphcitly 
that the 0{5) running couphng contributions are correctly resummed by the correction 
term in eq. (4.10): the leading-order term (4.13) in the expansion eq. (4.12) matches the 
naive (fixed coupling) dual eq. (4.7), while the first subleading correction matches the 
next-to-leading running coupling contribution eq. (4.8). 

We can thus write down the next-to-leading version of the resummed expression 
eq. (3.1): 

7(0;^, N) = 7/(q;s, N) + as^ii{as, N) + 0{al) - mom. sub.. (4-15) 
The leading-order term 77 is now given by 

7/(a„ N) = [a,7o(iV) + 7f (t) - ^] + 

+ 7$?)(«., N)-M, + [N - asc] - i/?o«. + ^ 0) - c 

(4.16) 

where 7^ is determined from eq. (4.4) using duality eq. (1.3), 7^^ is given by eq. (4.10), 
and the same values of the parameters Mg, k and c are used in these two contributions. 
Because, as discussed above, if k = kq and Mg — \ the resummed Xo eq. (4.4) reduces 
to the resummed xo of refs. [29,19], it follows that in such case ^^{ag/N) coincides with 
7s[q!s/(A'" — AA)] with AA — {c — coj/ag. In other words, this is the same as eq. (3.1), 
but with generic values of the parameters which characterize the behaviour of x about 
its minimum, eq. (3.3), and the cubic contributions to the Airy anomalous dimension 
accordingly included. 

The next-to-leading correction is 

7//(«a, N) = [asii{N) + 7f, (^) - eo - a. + 

1 /3oc /3o5/2\'/' N/as (4-17) 



— c 

a. 



Here, the term in square brackets is the analogue of the usual [19,20] next-to-lcading 
order term in the double-leading expansion (but determined from the new resummed x? 
eq. (4.4)) while the last two terms are the 0{a1/N) terms eq. (4.8) in the expansion of the 
Airy anomalous dimension. 

Let us now discuss the properties of this next-to-leading order expression. First, it is 
easy to see that the anomalous dimension eq. (4.15) only differs from the standard next-to- 
leading double-leading expansion of refs. [19-20] by sub-subleading terms. Indeed, we have 
seen that because of eq. (4.5) the reorganized expansion of the anomalous dimension in 
terms of 7^, 7^ . . . differs from the usual one by subleading terms. Furthermore, terms of 
order up to a1 in the expansion eq. (4.12) of the Airy anomalous dimension are explicitly 



17 



subtracted: note in particular that the sum of the order (3o term in 7/ and the order Pqc 
term in 7// exactly matches the subleading term eq. (3.12) in the expansion of ja- 

(3) 

The small behaviour of both 7/ and 7// is controlled by the pole of 7)^ at = 
A'^o, whose residue receives a 0{d) correction (compare eq. (4.10)), while singularities at 
N = asC in the expansion eq. (2.7) of the anomalous dimension are subtracted to the 
stated accuracy. As it is apparent from eq. (4.7), terms of increasingly higher order in 
the expansion of x about its minimum correspond to singular contributions to increasingly 
higher order derivatives of 7 at the point = asC. Namely, the quadratic term in x leads 
to a singularity in the first derivative of 7, the cubic and quartic term to a singularity in 
the second derivative, and so on. Hence, at the leading (7^) level, the square-root cut at 
N = OLgC is subtracted so that the anomalous dimension is continuous and has a continuous 
first derivative at N = ctgC, though it has discontinuous second derivative, which would be 
made continuous by including the quartic e term eq. (4.7) in the Airy resummation. 

At the next-to-leading (jss) level, thanks to the reorganization of the expansion 
eq. (4.4) there are no singularities in the fixed-coupling dual anomalous dimension nor 
in its first derivative. There are however singularities related to the running coupling term 
eq. (2.13): the pole (related to the quadratic contribution) and square-root singularity (re- 
lated to the cubic contribution) are subtracted, while the first derivative has a discontinuity 
that would be removed by including cubic corrections eq. (4.10) to 0{S^). 

At large N, the \/iV growth of 7^ is canceled by the 7^ subtraction, while the remain- 
ing constant large N limit is canceled by the 7^ subtraction, leaving a 1/VN drop. The 
linear and VN growth of 7^^^ are canceled by the 7^^^ and ^a\s subtractions respectively, 
leaving an 0{al) constant: 

„lL-(T<f'-7E-7!f,'„) = («.«^4- 

This constant may be included in the momentum subtraction term. This leaves then an 
0(1 /\/N) drop. Therefore, at large N the next-to-leading resummed anomalous dimension 
eq. (4.15) coincides with the standard GLAP anomalous dimension, since the subtracted 
small X eq. (2.7) terms drop at least as 1/A^ and the subtracted Airy terms drop at least 
as 1/\/N, while the GLAP anomalous dimension, as well known, grows (in modulus) as a 
power of In A". This property is preserved if the momentum subtraction is performed by 
subtracting the constant eq. (4.18), and then subtracting the remaining momentum non- 
conserving terms through eqs. (3.13)-(3.15). If momentum is instead enforced by simply 
subtracting a constant, the asymptotic large N resummed anomalous dimension will differ 
from the GLAP one by a fixed 0{a'^) constant. 

We conclude that to next-to-leading order the analytic properties of the leading order 
result eq. (3.1) are preserved: the branch cut at N = asC is subtracted from the small x 
anomalous dimension js+cus'Jss and replaced by the Airy pole, while the large N behaviour 
of the perturbative anomalous dimension 70 -|- as7i is preserved. 

We can now check that the next-to-leading order anomalous dimension eq. (4.15) 
differs from the leading order one eq. (3.1) displayed in figure 1 by a small correction. 
We expect this on the grounds that at large N they reduce to the corresponding GLAP 
anomalous dimension, while at small N for N ^ Nq they are controlled by the Airy 
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Figure 6: The improved anomalous dimension eq. (4.15) for UsC = 0.3 (lower dot- 
dashed), 0.4 (solid), 0.5 (upper dot-dashed), k = kq, Ms = 0.4 , g{N) eq. (3.15) with 
r = 1 and as = 0.2 compared with the NLO GLAP curve (dashed). The small kinks 
are due to cuts at N = asC which reappear at this accuracy through jss- 

pole, whose location depends very little on the input value of c. At next-to-leading order, 
however, we cannot use the perturbative values of the parameters c, k and Mg, because the 
next-to-leading order x does not have a minimum. This is due to the fact that, whereas 
the small M ~ drop of xi is removed by the double-leading resummation, fixed by 
momentum conservation, the large M ~ 1 drop of xi can only be removed through model- 
dependent assumptions, such as those of refs. [18,21]. Because the series of small M poles 
has alternating signs, this problem would disappear at the next-to-next-to-leading-log level, 
while at this order c, k and Mg must be treated as free parameters, upon the assumption 
that the minimum is restored by higher order corrections. 

The anomalous dimension eq. (4.15) is compared to its perturbative counterpart in 
figs. 6-8, for typical values of the parameters. It is seen that a reasonable stability of the 
results is obtained for 0.35 ^ Mg ^ 0.45, 0.3 ^ agC ^ 0.5, 3/4ko ^ k ^ 5/4ko, though at 
the extremes of the parameter range some instability due to large subleading corrections 
starts to appear. For instance, the discontinuity at = agC in the first derivative of the 
running coupling term eq. (2.13) starts showing up (figure 8) when n is lowered, since its 
coefficient is proportional to 

This shows that the leading order resummed expression eq. (3.1) is stable upon per- 
turbative corrections for a range of values of the parameters which are reasonably close to 
the leading order values eq. (3.4)-(3.6). However, me must conclude that the subleading 
correction eq. (4.17) cannot be used to improve the leading order resummation, because 
the values of the parameters which characterize it cannot be calculated in a reliable way. 
In fact, as already mentioned, the NLO perturbative X = Xo + <^sXi has no minimum. As 
a consequence, the uncertainty on this correction is of the same order as its size. 

Hence, the available information on the parameters c, k and Ms from the next-to- 
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Figure 7: The improved anomalous dimension eq. (4.15) for agC = 0.4, k = 5/4ko (lower 
dot-dashed), k = kq (solid), k = 3/4ko (upper dot-dashed), Ms = 0.4 , g{N) eq. (3.15) 
with r = 1, and as = 0.2 compared with the NLO GLAP curve (dashed). 




N 

Figure 8: The improved anomalous dimension eq. (4.15) for asC = 0.4, k = kq, Ms = 
0.35 (dot-dashed, lower at iV ~ 0.5), Ms = 0.4 (solid), Ms = 0.45 (upper dot-dashed), 
g{N) eq. (3.15) with r = 1, and Os = 0.2 compared with the NLO GLAP curve (dashed). 



leading order kernel xi is insufficient to improve the leading order resummed result. The 
only viable procedure, without further theoretical input, is to combine the leading-order 
resummation eq. (3.1) with the standard next-to-leading order correction to the GLAP 
equation, as we did in eq. (3.19) (figure 5). 
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5. Conclusion 

The main result of this work is the improved anomalous dimension given in eq. (3.1) 
(and its generalizations eqs. (3.18) and (3.19)). These relatively manageable expressions 
have the correct perturbative limit for small with fixed (which corresponds to large x 
for the splitting functions) and include the leading BFKL corrections with running coupling 
effects for small as with ag/N fixed. It is quite remarkable that the improved anomalous 
dimension in eq. (3.1), which does not contain free parameters, leads to a splitting function 
which is surprisingly close to the perturbative result down to small values of x (figure 2) in 
agreement with the data. The running coupling effects are essential to soften the small-x 
asymptotics, but a priori one would not expect that the corrected and the perturbative 
anomalous dimensions should be so close over an extended range at small x. The typical 
BFKL rise is softened and delayed to very small values of x^ while in the intermediate 
small-x region the splitting function is flat and close to the perturbative one. 

As is by now well known, a flat perturbative splitting function reproduces by evolution 
a rising structure function, and predicts a double logarithmic slope in quantitative agree- 
ment with that seen in the data [33]. However, the data suggest that the optimal value 
of the slope is contained in a region which is up to 15% smaller than the leading order 
GLAP value. It follows that the value of the splitting function in the flat region, being 
proportional to the square of the slope, can be at most 30% below its leading order GLAP 
value. Thus it is quite non trivial that our improved anomalous dimension directly fulfills 
this empirical constraint in the fiat region (where it lies below the perturbative splitting 
function) down to a; ~ 10~^. 

In comparison to the approach of refs. [18,21] we share the general physical framework, 
but there are two main differences. The first concerns the resummation procedure for 
x{q-s,M) near M = 0. We match systematically the perturbative expansion of x to that 
of the usual anomalous dimension, thereby finding a shift of order of the poles of x 
M = in a way which is made unambiguous by the constraint of momentum conservation. 
They instead, based on the change of scale variables in going from the BFKL symmetric 
into the DIS asymmetric configuration, displace the pole from M = to M + A^/2 = 0. 
In the relevant limit, ag small with ag/N fixed, the two displacements are of the same 
order, but in refs. [18,21] the input of momentum conservation, which controls our double- 
leading resummation, is not imposed although a posteriori the results are numerically 
compatible with it. Second, they determine the asymptotic small- a; behaviour by making 
some assumptions on the way to implement a symmetrisation M ^ (1 — M) of the BFKL 
kernel x(as,M), which we prefer not to do because of the model dependence it entails. 
For this we have to pay the price of some parameter dependence in the non leading terms. 
However, the fact that we find that after running coupling resummation the parameter 
dependence is weak gives support to their statement that the model dependence is not 
important. 

However, an advantage of our approach is that we end up with an explicit formal 
expression for the improved anomalous dimension which can be used in the conventional 
evolution equations and thus matched to standard perturbative evolution. In refs. [18,21] 
instead only the resummed BFKL kernel is determined analytically, so anomalous dimen- 
sions must be evaluated numerically from the solution. This makes the possibility of fitting 
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the data and matching to usual evolution equations more difficult. The recent work [34] by 
the same authors, which appeared while the present work was being completed, adopts a 
somewhat different symmetrisation procedure, such that, for example, momentum conser- 
vation is not automatic and is therefore imposed. A new result for the splitting function 
is found with a dip in the intermediate region at small x before the onset of the BFKL 
rise at even smaller values of x. This confirms the model dependence of the approach of 
refs. [18,21], which corresponds to the parameter dependence that we find. An advantage 
of our approach in this respect, however, is that the parameter dependence can be cleanly 
separated into the non-leading corrections. 

Indeed, the simple leading order expression of the improved anomalous dimension eq. 
(3.1) is only valid when non leading corrections in both the perturbative and the BFKL 
expansions can be neglected. We have discussed a general procedure for a systematic 
treatment of these higher order corrections, of which the NLO term is available, and we 
have considered it in detail. Following a general approach already used in our previous 
papers, our resummation strategy is to reabsorb in the leading term the parameters that 
describe the behaviour of xi^) ii^ar M = 1/2. In refs. [19,20] we resummed the value 
A = ttgC determining the small-x asymptotic behaviour as x~^. Here we include in the 
leading x(M) the first few terms of its expansion near the minimum at M = Mg that 
are needed for the Airy resummation. From the known non-leading correction we remove 
the perturbative expansion of the parameters included in the leading term. While the 
non leading correction is large before the shift we have shown that it reduces to a small 
correction after the subtraction. Thus we have shown that the non leading corrections are 
in principle under control. 

The difficulty is that the negative infinite behaviour of Xi(M) near M = 1 makes 
the treatment of the non leading contribution in the central region of M ambiguous. For 
example, the perturbative curvature k becomes negative in the central region when the 
Xi term is added, while we assume that the true curvature is positive (corresponding 
to a minimum). On the basis of the closeness of the leading-order result to standard 
perturbative evolution, and the phenomenological success of the latter, we expect that the 
subleading correction after subtraction is small. However, this also implies that it contains 
ambiguities of the same order as the correction itself. 

Thus, in practice, the known next-to-leading order correction is mainly useful as an 
estimate of the error on the leading approximation. The uncertainty on the non leading 
correction can be phcnomcnologically described in practice by replacing in the leading 
formula the lowest order perturbative values of the parameters describing the central region 
of X with generic values to be optimized to the data. For example, we have explicitly 
examined the variation of our results when the intercept c of the anomalous dimension 
is varied [eq. (3.18) and figure 3], and we also studied variations of the location of the 
minimum Mg and the curvature at the minimum k. 

In conclusion, we have achieved a precise quantitative understanding of the resum- 
mation of BFKL logarithms, which explains the empirical observation of the validity of 
GLAP evolution down to rather smaller values of x than might be naively expected. It 
will be interesting to compare parton distribution functions evolved with our resummed 
anomalous dimension with the latest HERA data to look for evidence [35] for the small 
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departures from GLAP evolution that we predict. 
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